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It is believed that soon after the Planck era, spacetime should have a semi-classical 

f^ ' nature. In this context we consider quantum fields propagating in a classical gravi- 

^D ■ tational field and study the backreaction of these fields, using the expected value of 

O . 

Q^. the energy-momentum tensor as source of the gravitational field. According to this 

^rjl theory, the escape from General Relativity theory is unavoidable. Two geometric 



counter-term are needed to regularize the divergences which come from the expected 
value. There is a parameter associated to each counter-term and in this work we 
found numerical solutions of this theory to particular initial conditions, for general 
Bianchi Type I spaces. We show that even though there are spurious solutions some 
of them can be considered physical. These physical solutions include de Sitter and 
Minkowski that are obtained asymptotically. 
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I. INTRODUCTION 

The semi-classical theory consider the backreaction of quantum fields in a classical geo- 
metric background. It began about forty years ago with the introduction of the quadratic 
terms in the Riemann tensor \^ made by the renormalization of the energy-momentum ten- 
sor, and continued with a more complete work made by De Witt J2|, and since then, its 
consequences and applications are still under research, see for example J3]. Some previous 
work 4] and 5] , studied the effect of particle creation in a Bianchi Type I spaces, and the 
renormalized coefficients of some of the higher derivative terms obtained by renromalization 
were set to 0. 

Different of the usual Einstein-Hilbert action, the predicted gravitational action allows 
differential equations with fourth order derivatives, which is called the full theory |1|, see 
also p. 

There are several problems connected to the full theory which prove that it is not con- 
sistent with our expectation of the present day physical world. Instability of flat space, 
Planck scale tidal forces, tachyonic propagation of gravitational particles, and violation of 
the positive energy theorem, see for example [if . Some of these results are obtained through 
linearization. 

n 

In the very interesting articles jtQ], Simon and Parker-Simon, suggest that the theory 
should be seen in a perturbative fashion. Their method results in differential equations for 
the metric to only second order. The order reduction allows for the elimination of the non 
physical effects. For instance in absence of a classical energy momentum source, Einstein's 
General Relativity is identically recovered in jq|. The above mentioned problems with the 
full theory are removed. The order reduction formalism should be the correct physical 
description in which the higher order terms are seen as a very small perturbative correction. 

also by Shapiro, Pelinson and others [1I|. Starobinsky idea was that the higher order terms 
could mimic a cosmological constant. In |lO] only the homogeneous and isotropic space time 
is studied. 

The full theory with four time derivatives is addressed, which apparently was flrst inves- 
tigated in Tomita's article [llj for general Bianchi I spaces. They found that the presence 
of anisotropy contributes to the formation of the singularity. Berkin's work shows that a 



quadratic Weyl theory is less stable than a quadratic Riemann scalar R'^. Barrow and Hervik 
found exact and analytic solutions for anisotropic quadratic gravity with A 7^ that do not 
approach a de Sitter space time. In that article Barrow and Hervik discuss Bianchi types // 
and Vlh and also a very interesting stability criterion concerning small anisotropics. There 
is also a recent article by Clifton and Barrow in which Kasner type solutions are addressed. 
H. J. Schmidt does a recent and very interesting review of higher order gravity theories in 



connection to cosmology 



id- 



In this present work, only the vacuum energy momentum classical source is considered 
in the higher derivative theory. It should be valid soon after the Planck era neglecting 
any particle creation that took place at that time. The full theory predicts explosions 
and formation of physical singularities depending in the parameters and initial conditions. 
From this point of view the full theory is certainly not a complete theory. In fact, the 
exact numerical solutions for Bianchi I spaces seem to reproduce all the problems mentioned 
above, this question was not investigated in this present work. 

But for some values of the parameters and very particular initial conditions, physically 
consistent solutions are obtained. The isotropization of spacetime also occurs with zero 
cosmo logical constant. The scalar Riemann four curvature oscillates near a constant value 
with decreasing amplitude. 

The following conventions and unit choice are taken R'^^^ = F"^^ — ..., Rab = Rad,^ -^ — ^a-i 
metric signature — h ++, Latin symbols run from — 3, Greek symbols run from 1 — 3 and 
G = h = c=l. 

II. FULL THEORY AND NUMERICAL SOLUTIONS 



The Lagrange function is. 



Metric variations in the above action results in 



k + R + a\ R^'Rab - ^^' ) + /5^' 



+ C,. (1) 



Eab = Gab +(/3- ^a) Hj^^ + aH^J - T,, - Ug^t, (2) 



where 

^ab = -^9abR — '^RRah + '^R-ah " '^'^Rdah, (3) 

^ab = -^SabRmnR"^"' + R;ab — '^R^^ Rcbna " D-Rafe " -^RQab-, (4) 

Gafe = Rab — -9abR, (5) 

and Tab is the energy momentum source, which comes from the classical part of the La- 
grangian Cc- Only classical vacuum solutions Tab = -^c = will be considered in this section 
since it seems the most natural condition soon after the Planck era. 

The covariant divergence of the above tensors are identically zero due their variational 
definition. The following Bianchi Type I line element is considered 

ds' = -dt' + e^'^^^'Ux^ + e^^^Wrfy^ + e'^'^^^'Uz^, (6) 

which is a general spatially flat and anisotropic space, with proper time t. With this line 
element all the tensors which enter the expressions are diagonal. The substitution of © in 
(j2)) with Tab = 0, results for the spatial part of (0), in differential equations of the type 

d^ . fd^ -A 

^a3 = /3^^a.,a.,d,J, (9) 

where the functions fi involve the derivatives of ai, 02, 03, in a polynomial fashion. The 
very interesting article |ll3| shows that the theory which follows from (^ has a well posed 
initial value problem. In llSll the differential equations for the metric are written in a form 

n 

suitable for the application of the theorem of Leray [1J| . 

Instead of going through the general construction given in [l3[, in this particular case, 
the existence and uniqueness of the solutions of Q-Q are more simply understood in |l5[. 

Besides the equations (Hjl-Q, we have the temporal part of (0). To understand the role 
of this equation we have first to study the covariant divergence of the equation ((21), 

VaE"^ = daE^^ + Tl^E'''' + r^,E"^ = 0, 

but since we are using ©-© to integrate the system numerically E'^^ = 0, then 

doE'' + r^o^°° + r[;oi^°° = 0. (10) 



Therefore if £"00 = initially, it will remain zero at any instant. Therefore the equation Eqq 
acts as a constraint on the initial conditions and we use it to test the accuracy of our results. 
For a space-like vector f " and a time-like vector t" = (1, 0, 0, 0) tidal forces are given by 
the geodesic deviation equations 

-^00/3 ~ ^otP (Oq, + da) ■ 

The theory predicted by (^ is believed to be correct if the tidal forces are less than 1 in 
Plank units units, 



l-Roaol — 1 (^° summation.^ 



(11) 



When this condition is not satisfied, quantum effects could introduce further modifications 
into dH). 

In order to integrate these equations we used a open source and well tested C library, the 
GSL GNU Scientific Library [l^, we used several algorithms provided by this library like 
Embedded Runge-Kutta-Fehlberg (4, 5) method and Implicit Bulirsch-Stoer method of Bader 
and Deuflhard, we also used Maple to calculate the equations and a Perl script developed 
by us to generate the C source code from the Maple output, and finally we used a AMD 
Athlon(TM) XP 2600+ to integrate the equations. 

A. Numerical Solutions 



For particular initial conditions and values for the parameters a, /?, A consistent numer- 
ical solutions for Q are shown in FIG. [Hand El In FIG. da = 0.1, (3 = -0.1, A = and 
in FIG. El a = 1000, /3 = -1000, A = 

The only non null initial conditions are chosen 



for FIG (Hand 
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FIG. 1: Using a = 0.1, (3 = -0.1, A = 
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FIG. 2: Using a = 1000, /3 = -1000, A = 
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FIG. 3: Using a = 0.1, f3 = 0.1, A = 
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FIG. 4: Using a = -0.1, /3 = 0.1, A = 
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FIG. 5: Using a = 0.1, /? = 0.1, A = 

for FIGI21 such that the 00 component of (j21) vanishes initially. According to pO|) . it is 
numerically checked that |-Eoo| < 10^^^ for the time interval in FICQand I-EqoI < 10^^^ for 
FIG. El 

The condition given in (fTT|) is satisfied initially, and during the time evolution in the 
solutions plotted in FIG. C]|21 Certainly these are physically accepted solutions since FIG. 
HI gives Minkowski space and FIG. |2 gives de Sitter space asymptotically. 

B. Explosions 



In FIG. lUthe values of the parameters are different 



a 



-0.1, /? = 0.1, A = 0. 



And again the only non null initial conditions are 
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are such that the 00 component of (J2I) vanishes. According to (jTUI). it is numerically checked 
that l^ool < 10^1^ for the time interval in FIG. El 
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It can be seen that although the initial conditions satisfy (fTTj). IRq^oI (^^ summation), 
assumes arbitrary large increasing positive values, which indicates the existence of explosions. 
In FIG. EJ it is shown that the curvature scalar increases to arbitrary large positive values. 

C. Formation of singularities 

The following values for 

a = 0.1, (3 = 0.1, A = 0, 

are the same for FIG. El and FIG. El In FIG.JHlthe non null initial conditions are 



and for FIG. E 
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such that the 00 component of (0) vanishes. According to ()10p. it is numerically checked that 
|ii^oo| < IQ-i^ for the time interval in FIG. Eland FIG. El Since both the scalar curvature 
tensor R and the squared Ricci tensor RabR'^'', increase abruptly, the solutions shown in 
FIG. Eland FIG. El are understood as singular type. The numerical error increases very fast 
when the solution comes closer to the singularity, which is expected. 

III. CONCLUSIONS 



In the present work it is considered general anisotropic Bianchi I homogeneous spacetimes. 
The full theory with four time derivatives problem are addressed. 

Only the vacuum energy momentum classical source is considered in the full theory. 
It should be valid soon after the Planck era and vacuum classical source seems the most 
natural condition. The full theory predicts explosions and formation of physical singularities 
depending in the parameters and initial conditions. From this point of view the full theory 
is certainly not a complete theory. In fact, the exact numerical solutions for Bianchi I spaces 
seem to reproduce all the problems mentioned above, this question was not investigated in 
this present work. 
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However, for some values of the parameters and very particular initial conditions, physi- 
cally consistent solutions are obtained. The isotropization of space time also occurs with zero 
cosmological constant. The scalar Riemann four curvature oscillates near constant values 
with decreasing amplitude. 

The above is understood as follows. The formation of singularities and Plank type explo- 
sions indicate that the theory certainly is not complete. The existence of physical consistent 
solutions show that the theory could have a space-time region of validity. 

Although we did not attempt a detailed verification, the numerical solutions obtained in 
this present work show no contradiction to the interesting previous calculations concerning 
anisotropics |ll|. The analytical solutions found by Barrow-Hervik and Clifton-Barrow can 
be understood as limit sets in the space of solutions of the quadratic theory (jT}. 
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